Prompt photon production and photon-hadron correlations at RHIC and the LHC 

from the Color Glass Condensate 
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We investigate inclusive prompt photon and semi-inclusive prompt photon-hadron production in 
high energy proton-nucleus collisions using the Color Glass Condensate (CGC) formalism which in- 
corporates non-linear dynamics of gluon saturation at small x via Balitsky-Kovchegov equation with 
running coupling. For inclusive prompt photon production, we rewrite the cross-section in terms 
of direct and fragmentation contributions and show that the direct photon (and isolated prompt 
photon) production is more sensitive to gluon saturation effects. We then analyze azimuthal cor- 
relations in photon-hadron production in high energy proton-nucleus collisions and obtain a strong 
suppression of the away-side peak in photon-hadron correlations at forward rapidities, similar to 
the observed mono-jet production in deuteron-gold collisions at forward rapidity at RHIC. We make 
predictions for the nuclear modification factor Rp(d)A and photon-hadron azimuthal correlations in 
proton(deuteron)-nucleus collisions at RHIC and the LHC at various rapidities. 

I. INTRODUCTION 



The Color Glass Condensate (CGC) formalism has been successfully applied to many processes in high energy 
collisions involving at least one hadron or nucleus in the initial state. Examples are structure functions (inclusive and 
^ [ diffractive) in Deeply Inelastic Scattering of electrons on protons or nuclei, and particle production in proton-proton, 
proton- nucleus and nucleus-nucleus collisions, for a recent review see Ref. The predicted suppression of RdA for 
the single inclusive hadron production in deuteron-nucleus (dA) collisions as well as the disappearance of the away 
O^l , side peak in di-hadron angular correlations in the forward rapidity region of RHIC 0, [1] are two of most robust 
^ ' predictions of the formalism which have been confirmed , see also Refs. 0-0| • The CGC formalism has been also 

successful in providing predictions for the first LHC data [g| in proton-proton (pp) and nucleus-nucleus collisions 
^y-^ ' [llj. Nevertheless, there are more recent, alternative phenomenological approaches which combine nuclear shadowing, 
transverse momentum broadening and cold matter energy loss to describe the RHIC data [13, Therefore, one 
needs to consider other observables which may help clarify the underlying dynamics of forward rapidity particle 
production at small x. 

Inclusive prompt photon production Ti"] and prompt photon-hadron angular correlations TF] in the forward rapidity 
region are two such examples. Furthermore, there are advantages to studying prompt photon production as compared 
to hadron production. It is theoretically cleaner; one avoids the difficulties involved with description of hadronization 
of final state quarks and gluons, usually described by fragmentation functions valid at high transverse momentum. 
Also, one does not have to worry about possible initial state-final state interference effects which may be present for 
5—1 I hadron production. In case of photon-hadron vs. di-hadron angular correlations, again the underlyin g th eoretical 
understanding is more robust. Unlike di-hadron correlations which involve higher number of Wilson lines jl6l |. photon- 
hadron correlations depend only on the dipole cross section properties of which are well understood. 

Both processes have been investigated previously, albeit not in detail and only in a limited kinematic range [l5| , 
see also Refs. pTl [l8|. In this work, we extend the existing results for inclusive prompt photon production by clearly 
separating the contribution of direct and fragmentation photons. We show that direct photons are more sensitive to 
gluon saturation effects in the kinematics regions considered. We then investigate the dependence of prompt photon- 
hadron azimuthal angular correlations on high gluon density effects and show that gluon saturation effects lead to 
disappearance of the away side peak. The effect is very similar to the disappearance of the away side peak in di-hadron 
correlations observed in the forward rapidity region of RHIC in dA collisions [l^. Therefore, a measurement of this 
correlation at RHIC and the LHC would greatly help to clarify the role of CGC in the dynamics of particle production 
at high energy. 

The advantage of the CGC formalism over the morephenomenological models is that the cross section for many 
of these processes have the same common ingredient [J [13, [Sj , the dipole total cross section; the imaginary part 
of the forward scattering amplitude of a quark-antiquark dipole on a proton or nucleus target. Its rapidity (energy) 
dependence is governed by the B-JIMWLK/BK evolutions equations [13, [13 and is pretty well- understood. The 
most recent advances in our understanding of the rapidity dependence of the dipole cross section include the running 
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coupling constant corrections and the full Next-to-Leading Order corrections [2J] . The only input is the dipole profile 
(dependence on the dipole size rt) at the initial rapidity which is modeled, usually motivated by the McLerran- 
Venugopalan (MV) model [1^. The sensitivity to this initial condition is expected to go away at very large rapidities, 
see Sec. Ill and Ref. ^. 

This paper is organized as follows; we consider prompt photon-hadron production cross section in section IIA and 
inclusive prompt photon production in section IIB where we describe how to separate the contribution of direct and 
fragmentation photons. We then present our detailed numerical results and predictions at kinematics appropriate for 
RHIC and the LHC experiments in section III. We summarize our results in section IV. 



II. THEORETICAL FRAMEWORK 



A. Semi-inclusive Photon-hadron production in proton-nucleus (pA) collisions 



The cross section for production of a quark and a prompt photon with 4-momenta I and k respectively (both 
on-shell) in scattering of a on-shell quark with 4-momentum p on a target (either proton or nucleus) in the CGC 
formalism has been calculated in [l^ and is given by 
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which, after using the explicit forms of the momenta in the expression for the trace, can be written as. 
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where the symbol T stands for a proton p or a nucleus A target, \fS is the nucleon-nucleon center of mass energj 



and Xq is the ratio of the incoming quark to nucleon energies such that p 



'3/2. The outgoing photon and 



quark rapidities are defined via k 



-^e^^ and I — ^e^' whereas A9 is angle between the final state quark and 

photon, cos{A9) = We note that this cross section was first computed in in coordinate space using the 

dipole formalism, the result of which agrees with the expression in Eq. Q after Fourier transforming to momentum 
space. 

The imaginary part of of (quark-antiquark) dipole-target forward scattering amplitude Np{bt,rt, Xg) satisfies the 
B-JIMWLK equation and has all the multiple scattering and small x evolution effects encoded. It is defined as 



NF{bt,rt,Xg) = ^ < Tr[l - V^{xt)V{yt)] >, 



(4) 



where is the number of color. The vector bt = {xl + yt)/2 is the impact parameter of the dipole from the target and 
ft = Xt ~ yt denotes the dipole transverse vector. The matrix V{yt) is a unitary matrix in fundamental representation 
of SU{Nc) containing the interactions of a quark and the colored glass condensate target. The dipole scattering 
probability depends on Bjorken Xg via the B-JIMWLK renormalization group equations. In the present case, it is 
related to the prompt photon and final state quark rapidities and transverse momenta via 
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(5) 



In order to relate the above partonic production cross-section to proton (deuteron)-target collisions, one needs 
to convolute the partonic cross-section in Eq. ([3]) with the quark and antiquark distribution functions of a proton 
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(deuteron) and the quark-hadron fragmentation function: 

-4- / dx,/K,g^) ; )> (6) 



d^foj dq^ dij-y drjhdO J zj^" ^/ J d^bt dk^ dtf drj-y drj^ dO 

where qt is the transverse momentum of the produced hadron, and f{xq,Q^) is the parton distribution function 
(PDF) of the incoming proton (deuteron) which depends on the hght-cone momentum fraction Xq and the hard scale 
Q. A summation over the quark and antiquark flavors in the above expression should be understood. The function 
Dh/q{zf,Q) is the quark-hadron fragmentation function (FF) where Zf is the ratio of energies of the produced hadron 
and quark ^ . Note that due to the assumption of coUinear fragmentation of a quark into a hadron, the angle A9 is 
now the angle between the produced photon and hadron. 

The light-cone momentum fraction related to the transverse momenta and rapidities of the produced 

hadron and prompt photon via (details are given in the appendix) 



1 A, „-r,. , <lt 



x„ = ^ kt e-"- + ^ 



= with = ^ (^-^ J . ,7) 



dHtdk^dijy (2^)3 kf J [zlt-ktY 



B. Single inclusive prompt photon production in proton- nuclear collisions 

The single inclusive prompt photon cross section can be readily obtained from Eq. ([T]) by integrating over the 
momenta of the final state quark. Integration over the quark energy is trivially done by using the delta function 
and leads to (after shifting It h — kt) , 

d^,i,)T^,iu)x ela.ra_^^^^_^^^__^^^^l_ I J ^^^ e^^''^--^NHbt,rt,x,), (8) 

where z = k^ /p^ denotes the fraction of the projectile quark energy p~ carried by the photon and drj^ = Various 
limits of this expression have been studied in |14| where it was shown that in the limit where photon has a large 
transverse momentum kt ^ zlt such that the collinear singularity is suppressed, one recovers the LO pQCD result 
for direct photon production process qg ^ qj convoluted with the unintegrated gluon distribution function of the 
target. On the other hand, if one performs the It integration above without any restriction, one recovers the LO pQCD 
expression for quark-photon fragmentation function convoluted with dipole scattering probability. In the limit where 
one can ignore multiple scattering of the quark on the target ("leading twist" kinematics), this expression reduces to 
the pQCD one describing LO production of fragmentation photons. It is therefore useful to explicitly separate the 
contribution of fragmentation photons from that of the direct photons. To this end, we rewrite Eq. (|S]) as 

-z 



d%td^ktdriy 7r(27r) 
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' ^^27rr^''^^ + i^-^)']]^ J d\ltNFixq,btJt), (9) 

where we have added and subtracted the second term. Notice that we use the same notation for coordinate rep- 
resentation of the forward dipole-target scattering amplitude Np and its two-dimensional Fourier transform. The 
second term in this expression describes production of direct photons whereas the first term gives the contribution of 



Since produced hadrons are assumed to be massless, we make no distinction between the rapidity of a quark and the hadron to which 
it fragments. Moreover, for massless hadrons, rapidity y and pseudo-rapidity rj is the same. 
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fragmentation photons. In order to see this more exphcitly, we let — > Zt + in the first term and keep the most 
divergent piece of the k integral to get 



(lu9(p)T^l(k) X ^^Fragmentation 



Direct 



cPbtcPktdrjj (PbtcPktdr]^ d^btd^kidr]^ 



= J^l^i/i^^^ fc* ) ^^^(^S' ^t' kt/z) + ^^^'[1 + (1 - ^f]^ dP^tl^t NF{xg,bt,lt) (10) 
where D^/q{z, fc|) is the leading order quark-photon fragmentation function (26j . 

D,/,{z,Q') = f|pi±ii^lnQVA^ (11) 

Eq. (jlOp is new and is our main result for single inclusive prompt photon production which includes contribution of 
both fragmentation (first term) and direct (second term) photons. In order to ensure that the divergence present in 
Eq. ([5]) is properly removed, one needs to regulate it self-consistently. Here we have done this separation by imposing a 
hard cutoff which would result in a mismatch between the finite corrections to our results and those that are included 
in parameterizations of photon fragmentation function, for example, using the MS scheme. However this mismatch 
is a higher order effect in the coupling constant and is therefore expected to be parametrically small. It should be 
noted that the separation between the direct and fragmentation contributions depends on the hard scale, chosen to 
be the photon transverse momentum, which is already well-known in pQCD. 

Eq. (fTO)) exhibits some interesting features; the dipole scattering probability Np is probed at kt/z (where kt is 
the external momentum) in case of fragmentation photons whereas it depends on the internal momentum It in case 
of direct photons. Furthermore, in case of direct photons, the integrand is peaked at values of transverse momenta 
h ~ Qs- This means that fragmentation photons should be much less sensitive to high gluon density effects than 
direct photons since they probe the target structure at higher transverse momenta. This will be verified numerically 
in the following sections. 

In order to relate the partonic cross-section given by Eq. (|10p to photon production in deuteron (proton)-nucleus 
collisions, we convolute Eq. ([TO]) with quark and antiquark distribution functions of the projectile deuteron (or proton), 

dxMx,,k^) + f,{x,,kf)]^^-^—^^ , (12) 



d'^btd'^ktdrjj J^'^in d'^btd'^ktdij.y 

where a summation over different fiavors is understood. Equations (|10ll2p are our final results for the single inclusive 
prompt photon production. The light-cone fraction variables Eq. (|10ll2p are defined as follows. 



2 ^ g 
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I,- i- mm . 

^ ^ fc_ ^ = ^ ^ith = z™„ = -^e^^. (15) 

P XqVS Xq y' S 

where in Eq. (|14p the right hand-side approximation is valid if It ^ kt- Notice that since now Xg depends on the angle 
between It and kt , the integral over the angle in Eq. ([TU)) is not more trivial and can be done numerically. One should 
also note that the light-cone fraction variables defined above for the inclusive prompt photon cross-section Eqs. (|10|12p 
are different from the corresponding semi-inclusive hadron-photon cross-section Eqs. (|3l6p defined in Eqs. ([7]), see the 
appendix for the derivation. 



III. NUMERICAL RESULTS AND PREDICTIONS 



The forward dipole-target scattering amplitude appears in both semi-inclusive photon-hadron and inclusive prompt 
photon cross-section Eq. ()6I12I) and incorporates small-x dynamics which can be computed via first principle non- 
linear B-JIMWLK equations |22j in the CGC formalism. In the large Nc limit, the coupled B-JIMWLK equations are 
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simplified to the Balitsky-Kovchegov (BK) equation \2^, a closed-form equation for the evolution of the dipole am- 
plitude in which both linear radiative processes and non-linear recombination effects are systematically incorporated. 
The running-coupling BK (rcBK) equation has the following simple form: 

' ' ^ d^n K'''''{r,ri,r2)[NF{ri,x) + NF{r2,x) - NF{r,x) - NF{ri,x)NF{r2,x)] , (16) 



d\n{xo/x) 

where r2 = f—fi. The evolution kernel X''"" is given by Balitsky's prescription [27] with the running coupling. 
The explicit form of iC"" with details can be found in Refs. [13, [lal- The only external input for the rcBK non- 
linear equation is the initial condition for the evolution which is taken to have the following form motivated by 
McLerran-Venugopalan (MV) model [25| . 



Af{r,Y = 0) = 1-exp 



Ar 



(17) 



where A = 0.241 GeV [IS H!]. The initial saturation scale Qos (with s — p^A for a proton and nuclear target) at 
starting point of evolution (at Xq = 0.01) and the parameter 7, are free parameters which are determined from a fit to 
other experimental measurements at small-x. It was shown that inclusive single hadron data in pp collisions at RHIC 
can be described with a initial saturation scale within Q^^ = 0.168 0.336 GeV^ [1, H^]. However, HERA data 
on proton structure functions prefers the lower value for the proton initial saturation scale Qg^ « 0.168 GeV^ [29| . 
One have also freedom to run 7 as a free parameter in the minimization and obtain its preferred value in a fit to 
HERA data. In order to investigate the uncertainties due to initial condition of the rcBK equation, we will consider 
the following three parameter sets which all provide an excellent fit to the HERA data for proton targets [13, : 





set I : 




= 0.2GeV^ 7 = 1, 




set II : 




= 0.168 GeV^ 7 ==1.119, 




set III : 




= 0.157 GeV^ 7 = 1.101. 


In the MV model [2 


5], the parameter 7 in 


Eq. ini) 


is 7 = 1. However, it has been recently shown 



(18) 

[m that the 

effective value of 7 can be larger than one when the sub-leading corrections to the MV model are included [l^]. The 
parameter 7 appears to be also important in order to correctly reproduce the single inclusive particle spectra, and a 
larger value 7 > 1 is apparently preferable at large-fct (gI. [lOj. 

In our approach the difference between proton and nuclei originates from different initial saturation scales Qos in the 
rcBK equation via Eq. I|17p. In the case of inclusive hadron production in proton- nucleus collisions, due to theoretical 
uncertainties and rather large errors of the experimental data, it is not possible to uniquely fix the initial value of 
Qqa- In the case of minimum- bias dAu collisions, the initial nuclear (gold) saturation scale within Qq^ = 3-^4 Q^^ is 
consistent with the RHIC inclusive hadronproduction data [3, [E S HO] ■ The extracted value of Qqa is also consistent 
with the DIS data for nuclear targets [1^, [sO]. Here, we will also consider the uncertainties due to the initial condition 
of the rcBK equation for a nuclear target. Note that Qqa should be considered as an impact-parameter averaged value 
since it was extracted from the minimum-bias data. For the minimum-bias collisions, one may assume that the initial 
saturation scale of a nuclei with atomic mass number A, scales linearly with A^/'^ , namely we have Qg^ = cA^/'^ Qq^ 
where the parameter c is fixed from a fit to data. In Ref. |30|], it was shown that NMC data can be described with 
CSS 0.5. 

We wiU use the NLO MSTW 2008 PDFs [H and the NLO KKP FFs [H]. For the photon fragmentation function, 
we will use the full leading log parametrization [2^, [s^ . We assume the factorization scale Q in the FFs and the 
PDFs to be equal and its value is taken to be qt and kt for the semi-inclusive and inclusive prompt photon production, 
respectively. 



A. Direct and fragmentation prompt photon in pp and pA collisions at RHIC and the LHC 

We start by considering direct and fragmentation photon production in pA collisions at RHIC and the LHC. In 
nuclear collisions, nuclear effects on single particle production are usually evaluated in terms of ratios of particle 
yields in pA and pp collisions (scaled with a proper normalization), the so-called nuclear modification factor RpA- 
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FIG. 1: Nuclear modification factor for direct plioton production in minimum-bias pA (daslied line) and dA (solid line) collisions 
at RHIC {VS = 0.2 TeV) at = 2, 3. The curves are obtained from Eq. (fT2)l using the solution to rcBK equation with the 

?0p 



initial saturation scale Qq^ — 0.168 GeV^ for proton and Qqa — 3Qop for £^ nucleus (gold). 



The nuclear modification factor Rp(d)A is defined as 



_ 1 dN'^^^'^^ dNPP^'^^ 
''^ " 2N,oU cPpTdrj ' d^pTdv ' 
_ 1 dNP^^-f^ dNPP^-^^ 

where the yield '^^''^^^p^^'^ ^ can be calculated from the invariant cross-section given in Eq. (IT^ . The normalization 
constant Ncoii is the number of binary proton-nucleus collisions. We take Ncou = 3.6, 6.5 and 7.4 at = 0.2, 4.4 and 
8.8 TeV, respectively [S^. In order to compare our predictions for it!^^ with the experimental value, one should take 
into account possible discrepancy between our assumed normalization Ncoii and the experimentally measured value 
for Ncoii by rescaling our curves. Again we expect that some of the theoretical uncertainties, such as sensitivity to K 
factors (which effectively incorporates the missing higher order corrections), will drop out in .Rp^^j^^. 

In Fig. [l]we show the nuclear modification factor for both pA and dA collisions at RHIC. This is to facilitate a 
comparison of and to distinguish between the genuine saturation effects in the target nucleus from isospin effects 
in the projectile deuteron. Clearly there is a large difference between a proton and a deuteron projectile as far as 
prompt photon production is concerned. This difference is more pronounced in the forward rapidity region and at 
high transverse momentum where one probes the quark content of the projectile. This is due to difference between 
the up and down quark distributions of a proton (note that nuclear effects in the wave function of a deuteron 
are ignored as they are known to be small). This is a well known effect, nevertheless, for sake of clarity and to 
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FIG. 2: Nuclear modification factor for direct, fragmentation and inclusive prompt photon production in minimum-bias p(d)A 
collisions at RHIC VS = 0.2 TeV (right) and the LHC VS = 4.4 TeV (left) energy at various rapidities. The curves are the 
results obtained from Eq. (|12|l and the solution to rcBK equation with the initial saturation scale Q^p = 0.168 GeV^ for a 



proton and Qqa = 3Qop for a nucleus (gold), corresponding to set II in Eq. ([18] 



illustrate the difference between inclusive hadron production and QED probe namely prompt photon production, 
we illustrate this in some detail. In case of photon production, the production cross section is weighed by the the 
charged squared of a given quark flavor. For example (assuming only two flavor), for a proton projectile this is given 
by /q/p — (2/3)^ Up + (1/3)^ dp where Up, dp denote the distribution functions of up and down quarks in a proton. 
Ignoring nuclear effects in a deuteron, we assume a deuteron is a system of free proton and a neutron in which case the 
corresponding expression is e'^ fqf^ = (2/3)^ + (1/3)^ dp + (2/3)^ m„ + (1/3)^ d„ where u„,d„ denote the distribution 
functions of up and down quarks in a neutron. Assuming isospin symmetry gives u„ = dp and dn — Up which leads to 
(5/9) [up + dp] for a deuteron. Comparing this expression with two times that of a proton, the relative contribution 
of up quarks in a deuteron (5/9) is smaller than that of a twice a proton (8/9). Since there are more up quarks than 
down quarks (by a factor of 2 3 in this kinematics) in a proton, and their relative weight is smaller, this leads to a 
further reduction of RdA as compared with RpA in prompt photon production, see Fig. [TJ We note that in the absence 
of the charged squared factor, which is the case for inclusive hadron production, one would get d = p + n = 2p as one 
should since possible nuclear effects in a deuteron are ignored here. At the LHC the isospin effect is absent since the 
same projectile is used for the pp and pA collisions. This helps to understand the physics of QCD saturation more 
clearly, as the suppression of the signal will not be contaminated with isospin effect. 

In Fig. [5J we show the minimum-bias nuclear modiflcation factor for the direct, fragmentation and the inclusive 
prompt photon production at RHIC and the LHC energies ^/S = 0.2,4.4 TeV at various rapidities 77 obtained from 
Eqs. ([T2I [191) supplemented with rcBK solution Eq. ([T6|) with the initial saturation scale for proton Qq^ « 0.168 GeV^ 
and nuclei Q04 = SQop- It is seen that the nuclear modiflcation Rjf^^-^j^ for the fragmentation photon is bigger than 
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FIG. 3: Nuclear modification factor for direct photon production in p(d)A collisions at various rapidities at RHIC y/S = 0.2 
TeV (right) and the LHC y/S — 4.4 TeV energy (left). The curves are the results obtained from Eq. p2p and the solution to 
rcBK equation using different initial saturation scales for a proton Qop and a nucleus Qoa- The band shows our theoretical 
uncertainties arising from allowing a variation of the initial saturation scale of the nucleus in a range consistent with previous 
studies of DIS structure functions as well as particle production in minimum-bias pp, pA and AA collisions in the CGC 
formalism, see the text for the details. 



the direct and inclusive prompt photon. This is what we expected in our picture since direct photon cross-section in 
Eq. (flO)) probes the target structure function at lower transverse momentum kt (and consequently lower x) than the 
fragmentation part with transverse momentum kt/z and therefore is more sensitive to the suppression of structure 
function and the saturation effect. However, as we increase the energy the enhancement of the fragmentation photon 
-^'p{d)A RHIC will be also replaced with suppression at the LHC, see Fig. [2] top panel. This is simply due to the 
fact that both the fragmentation and the direct part Eq. (jlOp depend on the color dipole forward amplitude which 
encodes the small-x dynamics and at higher energy, the small-x evolution leads to suppression of Rj^^^^^j^- 

In a collider experiment such as the LHC, the secondary photons coming from the decays of hadrons, overwhelm the 
inclusive prompt photon measurements with order of magnitudes. In order to reject the background, isolation cuts 
arc imposed [33|- Contribution of fragmentation prompt photon is reduced by imposing an isolation cut^. A proper 



^ If we assume that pc is the total transverse momentum of a fragmentation jet, the photon's energy is then Ej = zpc and the total 
hadronic energy within the jet is Ef^ = (1 — z)pc. By isolation cut criterion, the hadronic energy does not have to be more than eEj in 
the isolation cone. This gives the lower limit of z (or Xq convolution) in Eq. II13I I integral, namely Zc = + t) < z. Given that the 
integrand of fragmented part is proportional to I/2 and dominated at lower limit of integrand, we expect that the isolation cut reduces 
the fragmentation contribution more severely than the direct one. 
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FIG. 4: Nuclear modification factor for direct photon (right) and inclusive prompt photon (left) production in pA collisions 
at various rapidities a the LHC y/S — 8.8 TeV energy. The band (CGC-rcBK-av) similar to Fig. [3] corresponds to the results 
obtained from Eq. (|12l) and the solutions to the rcBK evolution equation using different initial saturation scales for a proton 
Qop and a nucleus Qoa, see the text for the details. 



incorporation of the isolation cut criterion in our framework is beyond the scope of this paper. However, from Fig. [2] 
it is seen that at higher energy at forward collisions, R'^f^^-^j^ for direct and single inclusive prompt photon becomes 
remarkably similar, indicating that to a good approximation, one may assume that the nuclear modification factor 
for direct and isolated prompt photon are equal. 

In Fig. 131 we show the minimum-bias nuclear modification factor for the direct photon production at RHIC and the 
LHC energies yfS = 0.2, 4.4 TeV at various rapidities rj obtained from rcBK solutions Eq. ([T5)) with the initial proton 
saturation scale Qop ~ 0.168 and 0.2 GeV^ corresponding to parameter sets I and II in Eq. (fT8)) . For nuclear target 
in minimum-bias collisions, we take two initial saturation scales for nuclei ( gold and lead) = 3 ^ ^Q^p which are 
extracted from a fit to other experimental data on heavy nuclear target [sl. la. l30l| . For a proton target, we have checked 
that parameter sets II and III give similar results for Rp(^^^j^ with better than 10% accuracy. Therefore, in Fig. |3]we 
only show results obtained from two parameter sets I and II in Eq. psp . The band in Fig. |3] shows our uncertainties 
arising from a variation of the initial saturation scale of the nucleus in a range consistent with previous studies of DIS 
structure functions as well as particle production in minimum-bias pp, pA and AA collisions in the CGC formalism. 
One may therefore expect that the possible effects of fluctuations (of nuclcons in a nucleus) on particle production is 
effectively contained in our error band. 

From Fig. [3l it is seen that the nuclear modification for direct photon production is very sensitive to the initial 
saturation scale in proton and nuclei. However, this uncertainties will be reduced for more forward collisions at 
higher energy at the LHC. The same effect has been observed for the inclusive hadron production in pA collisions @]. 
This clearly indicates that the nuclear modification in p(d)A collisions is a sensitive probe of saturation effects and 
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FIG. 5: Right: Nuclear modification factor for direct photon production at ?7 = 3 in minimum- bias dAy/S = 0.2 TeV (RHIC) 
and pA = 4.4, 8.8 TeV (LHC) colUsions. The curves are the results obtained from Eq. (|12p and the solution to rcBK 
equation with the initial saturation scale Qq^ — 0.168 GeV^ for a proton and Qqa ~ 3Qop ^or a nucleus. Left: Comparison of 
the inclusive prompt photon nuclear modification factor predictions from the CGC (in this paper) and the standard coUinear 
factorization approach p^ . The band CGC-rcBK-av is the same as in Fig. 2] 



^p{d)A measurements for direct photon and inclusive hadron provide crucial complementary information about initial 
saturation scale and small-x evolution dynamics. In Fig. [31 it is seen that at a fixed rapidity and energy for a fixed 
initial saturation scale for proton Qop, a bigger initial saturation scale for nuclei Qqa leads to a bigger broadening 
and consequently enhances the cross-section and -Rp(d)^ if Ncou is kept fixed. 

In Fig. m we show our predictions for RpA for direct photon (right) and inclusive prompt photon (left) production 
in pA collisions at various rapidities a the LHC VS = 8.8 TeV energy. The band (CGC-rcBK-av) similar to Fig. [3] 
corresponds to the results obtained from Eq. (fT2|) with the solutions of the rcBK evolution equation ([T6l) . 

In Fig. [5] (right), we compare Rpf^j^j^ for direct photon at 77 = 3 for RHIC energy ^/S = 0.2 TeV and the LHC 
energies 4.4, 8.8 TeV. It is seen that the suppression of R^^^ at the LHC is larger compared to RdA at RHIC and 
persists at higher transverse momentum. This larger suppression is even more impressive given that fact that a good 
amount of the observed suppression of RdA at RHIC is due to the projectile being a deuteron rather than a proton. 
In Fig. [S] (left), we compare the CGC prediction (CGC-rcBK-av) obtained here with the collinear factorization result 
(EPS09) [I4I for inclusive prompt photon Rjj^ at 77 = 3 at the LHC. It is seen that the LHC measurements of the 
inclusive prompt photon at forward rapidities can discriminate between the collinear (standard parton model) and 
the CGC approach. 

Some words of caution are in order here. Strictly speaking our formalism is less reliable for collisions at around 
mid-rapidities and high transverse momenta. This is due to the fact that our formula is valid for asymmetric collisions 
like pA or pp collisions at forward rapidities when a projectile can be treated in the standard collinear approximation 
while for the target we systematically incorporated the small-x re-summation (at the leading log approximation) 
effects. Note however, for the case of pp collisions (our reference for RJa) RHIC, the saturation scale of target 
proton is rather small, and it is not clear that the CGC formulation will be applicable. Moreover, our parameter sets 
in Eq. ((THl) was obtained from a fit to HERA data at small-x x < 0.01 and for virtualities G [0.25, 40] GeV^ [29j . 
Therefore, our predictions are less reliable at high-fct {kt > 6 7 GeV). 
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FIG. 6: The relative azimuthal correlation P{A9) defined in Eq. ((20} for minimum-bias p(d)A and pp collisions at RHIC 
y/S — 0.2 TeV (upper) and the LHC \/S = 4.4 TeV energy (lower) obtained from the rcBK solutions with different initial 
saturation scales. 



B. Prompt photon-hadron correlations at RHIC and the LHC; the signature of saturation 



We now focus on azimuthal angle AO correlations of the prompt photon-hadron spectrum, where the angle A9 
is the difference between the azimuthal angle of the measured hadron and single prompt photon. We present our 
predictions for semi-inclusive prompt photon-hadron (for hadron we consider only neutral pion here) production at 
RHIC and the LHC in pp and p(d)A collisions in terms of P{A6) defined as follows, 

J p(<i) T-^h{q) 7(fe) X J p{d) T-^h{q) 7(fc) X 

P{M)^ — [AO] I — {AQ = A6'J , (20) 

<Pht dkf dq^ dy-y dyi dO d^ht dkf dqf dy-y dyi dO 

where the prompt photon-hadron cross-section in above expression is given in Eq. This definition has a simple 
meaning of the probability of, the single semi-inclusive prompt photon-hadron production at a certain kinematics and 
angle AO given the production with the same kinematics at a fixed reference angle AOc- We take AOc = 7r/2. As 
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FIG. 7: The relative azimuthal correlation P{A9) defined in Eq. (f20|l for minimum-bias dA collisions at RHIC VS = 0.2 TeV 
for two different windows of kinematics of transverse momentum of produced prompt photon kt and hadron (neutral pion) qt 
at fixed rapidity rih = r/j = 3. The curves are the results obtained from the rcBK equation solution with the initial saturation 



scale Qop = 0.168 GeV for proton and Qq 



3<3op for a nuclei. 



will show the P{A6) defined in this way has a non-trivial structure and can probe the physics of small-x and gluon 
saturation. In principle, one is free to chose a different reference angle A.9c, however any value A9c << tt will only 
change the normalization rather than the main picture. The advantage of the above definition for the azimuthal 
correlations is that it is experimentally easier to measure as it does not require a different experimental setup and 
run for the trigger or reference. Moreover, in dA collisions at RHIC, the isospin effect in P{A9) will drop out via 
normalization'^ and this facilitates to single out the importance of the saturation effect at forward rapidities in contrast 
to the nuclear modification factor -R^^. 

In this approach a fast valence quark from the projectile proton radiates a photon before and after multiply 
scattering on the color glass condensate target, see Fig. [SI In this picture, the projectile is treated in the collinear 
factorization, and therefore the photon radiation from quark at this level has the standard features of pQCD, including 
the back-to-back correlation in the transverse momentum. As a result of multiple scatterings, the quark acquires a 
transverse momentum comparable with the saturation scale, the only relevant scale in the system, and the intrinsic 
angular correlations are washed way. 

In Fig. ini we show P{A9) at forward rapidity 77/1 = 77^ = 3 for qh = fcf = 2 GeV at RHIC and qh — kt — 6 GeV at 



^ We checked that numerically the isospin effect brings less than 2% contribution to the azimuthal correlation defined via Eq. II20I I. 
Therefore, due to our particular definition of P{A9) in Eq. I I20II . the differences between a deuteron and a proton projectile are negligible 
unlike the prompt photon production case. 
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FIG. 8: The relative azimuthal correlation P{A9) for minimum-bias p(d)A collisions at different energies VS = 0.2,4.4, 8.8 
TeV for a fixed rapidity rih = ri-y = 3 and transverse momentum of the produced prompt photon kt = 5 GeV and hadron 
(neutral pion) qt — 3 GeV. The curves are the results obtained from the rcBK equation solution with the initial saturation 
scale Qop = 0.168 GeV^ for proton and Qoa ~ 3Qop for a nuclei 



the LHC for two different initial saturation scale for proton Qq^ = 0.168, 0.2 GeV^ and nuclei Qq^ = 3 ^Qtp- For 
such low ptS we are most likely probing the saturation region of the nuclear wave function due to the small values 
of Xg. It is clear that the away-side prompt photon- hadron cross-section (at « tt) is suppressed for the bigger 
saturation scale (corresponding to a denser system). It is also seen from Fig. [5] that P(A6') is very sensitive to the 
initial saturation scale. Unfortunately, this bring rather large theoretical uncertainties. However, as we will show 
in the following the suppression of the way-side correlations seems to be a robust feature of our results and it less 
depends on our theoretical uncertainties. 

In Fig. [7] we show the relative azimuthal correlations obtained from the rcBK solution for a fixed initial saturation 
scale = 0.168 GeV^ and Qg^ = '■^Qop at forward rapidities rjh = rj-y — aX RHIC ^/S = 0.2 TeV for two different 
kinematics windows of transverse momenta: We show in top panel, the results with a fix transverse momentum of 
prompt photon qt = 5 GeV at different transverse momentum of produced hadron, and in down panel, with a fixed 
transverse momentum of hadron fc^ = 5 GeV but at various transverse momentum of the produced prompt photon. 
It is seen, the relative azimuthal correlation is suppressed at A9 = tt as the transverse momentum of the produced 
hadron or prompt photon decreases and becomes comparable to the actual saturation scale of the system. This is the 
case regardless which of two transverse momenta of the hadron or prompt photon decreases. In Fig. [7] lower panel, it is 
seen that when the prompt photon transverse momentum becomes comparable with the saturation effect the away-side 
azimuthal angular correlation of photon-hadron completely washes away. The same effect happens at lower transverse 
momentum of the produced hadron. This is simply because of fragmentation effect namely the transverse momentum 
of the produced parton (that should be compared with the saturation scale) is higher than the transverse momentum 
of the fragmented hadron. Again, the suppression of away-side correlations is clearly due to the saturation effect since 
as we lower the transverse momentum of the produced particle, the system of hadron-photon become more sensitive to 
the small-x gluon saturation. Note that the hadron-photon cross-section in Eq. ([6]) has collinear singularity. Therefore 
for a proper investigation of the correlations at sa 0, in principle, one should first extract the collinear singularity 
in a same fashion as demonstrated in section II by introducing the quark-photon fragmentation function. Therefore, 
our results at near-side AO « should be less reliable. Nevertheless, we expect that the sensitivity to the collinear 
singularity effect should drop out in the correlation defined in Eq. (|20|) via normalization. We checked that contrary 
to the away-side correlations, the near-side peak is not sensitive to the saturation physics as the correlations does not 
change with varying the density of the system, see also Fig. [SI 

In order to further understand the relative sensitivity of the away side peak to saturation dynamics, in Fig. [51 
P{A9) we compare at various energies at RHIC and the LHC for a fixed transverse momentum of the produced 
prompt photon kt = 5 GeV and hadron qt — 3 GeV at rapidity rjh — rj^ = 3. It is clear that the away side peak 
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goes away as one increases the energy. Again this is due to the fact that as we increase the energy, the gluon density 
increases and non-hnear gluon recombination or the saturation effect becomes important. From Fig. [SJ it is obvious 
that at the LHC the away-side azimuthal correlations of photon-hadron will be strongly suppressed. 

We conclude that the suppression of the away-side azimuthal photon-hadron correlations defined via Eq. (|20p . 
with decreasing the transverse momentum of the produced prompt photon or hadron, or increasing the energy, or 
increasing the size/density of system, all uniquely can be explained within the universal picture of gluon saturation 
without invoking any new parameters or ingredients to our model. 



IV. SUMMARY 



We have investigated prompt photon production and prompt photon-hadron azimuthal angular correlations in 
proton-proton and proton-nucleus collisions using the Color Glass Condensate formalism. We have provided predic- 
tions in the kinematic regions appropriate to RHIC and the LHC experiments. We have shown that single inclusive 
and direct prompt photon production cross section in p(d)A collisions at forward rapidities at both RHIC and the 
LHC is suppressed, as compared to normalized production cross section in proton-proton collisions. At RHIC, a good 
portion of the predicted suppression is due to the projectile being a deuteron rather than a proton. This suppression is 
larger at the LHC compared to RHIC which is even more impressive given that the projectile at the LHC is a proton. 
We showed that direct photon production is most affected by gluon saturation effects in the target nucleus than the 
fragmentation photons. However, at the LHC energies at forward rapidities the nuclear modification suppression for 
direct, fragmentation and inclusive prompt photon production is rather similar. We showed that the nuclear modi- 
fication factor fip(^)yj for inclusive prompt photon production at RHIC and the LHC is a sensitive probe of small-x 
dynamics. We note that our results based on gluon saturation dynamics and using the Color Glass Condensate for- 



malism are different from those coming from the collinear factorization [12|. Therefore, Rj^^^j^ measurement at RHIC 
and the LHC is a crucial test of different factorization schemes, see also Refs. IH, for other observables. 

We have also studied prompt photon-hadron azimuthal angular correlations in kinematic regions which can be 
probed by RHIC and the LHC experiments. It is shown that the away side peak in photon-hadron angular correlation 
goes away as one lowers the final state particle's momenta, very similarly to the disappearance of the away side peak 
in di-hadron correlations in forward rapidity dA collisions at RHIC. At fixed transverse momenta, the suppression of 
the away side peak gets stronger as one goes to larger rapidities (more forward) or higher energy or denser system 
as expected, due to stronger saturation effects in the target nucleus. Presently, wc are not aware of any alternative 
approach which leads to this novel phenomenon. Note that in contrast to the nuclear modification factor for prompt 
photon, the prompt photon-hadron azimuthal angular correlation defined via Eq. (1201) is free from the isospin effect, 
and can be considered as a cleaner probe of saturation effect. Finally, we emphasize that prompt photon-hadron 
azimuthal angular correlations suffers from much less theoretical uncertainties as compared to di-hadron azimuthal 
angular correlations and thus a measurement of this correlation would go a long way toward establishing the dominance 
of gluon saturation effects at small Xg. 

It will be interesting to see what the predictions of pQCD-motivated models [l3| are for photon-hadron azimuthal 
angular correlations. In these models one usually needs to combine models of higher twist shadowing, the Cronin 
effect and cold matter energy loss in order to describe the data on single inclusive hadron production and di-hadron 
azimuthal angular correlations. The advantage of the CGC formalism is that the same framework can be used to 
describe nuclear shadowing of structure functions 40] at small x and includes transverse momentum broadening (the 
Cronin effect) [4l| . It does not however include cold matter energy loss due to longitudinal momentum transfer 
between the projectile and the target which may be important at the very forward rapidities. It is not clear at the 
moment how to calculate this effect from first principles QCD. Even though this energy loss itself is small, due to 
steepness of the production cross section at forward rapidity, it can suppress the cross section significantly. 



Appendix A 

The purpose of this appendix is to define the kinematics and derive the needed relations between various light-cone 
energy fractions which appear in the production cross sections used. This is slightly different from the standard 
relations used in production cross sections based on collinear factorization theorems of pQCD. We first consider 
scattering of a quark on the target where a photon and a quark are produced, depicted in Fig. 1, 

q{p)A{pA)-^l{k)q{l)X, (Al) 

where A is a label for the multi-gluon state, described by a classical field representing a proton or nucleus target. In 
the standard pQCD (leading twist) kinematics, only one parton from the target interacts. This is not the case here 
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FIG. 9: The diagrams (at leading log approximation) contributing to the prompt photon+hadron production within the color 
glass condensate picture. The crossed white blob denotes the interaction of the projectile quark to all orders with the strong 
background field of the target nucleus A. The black blob represent the quark-hadron fragmentation process. 



since the target is described by a classical gluon field representing a multi-gluon state with intrinsic momentum rather 
than an individual gluon with a well defined energy fraction Xg and zero transverse momentum. Nevertheless, since 
most of the gluons in the target wave function have momentum of order Qs, one can think of the state describing the 
target as being labeled by a (four) momentum pA- In this sense, the gluons in the target collectively carry fraction 
Xg of the target energy and have intrinsic transverse momentum denoted PA.t- This also means that there is no 
integration over Xg in our case unlike the coUinearly factorized cross sections in pQCD (this basically corresponds to 
setting Xg equal to the lower limit of Xg integration in pQCD cross sections). We thus have 
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qf^ = 






(k-, 




1^/21-, h) , 



r 



(A2) 



where , P^, arc the momenta of the incoming projectile, target and the produced hadron respectively. (Pseudo)- 
rapidities of the produced quark and photon are related to their energies via 



(A3) 

- k) and using Eq. (|A2I) leads to 

(A4) 
(A5) 
(A6) 

The above relations and Eq. (|A2I) (and the on mass shell condition) can be used to derive the following expressions 
for the energy fractions Xq,Xg. We obtain, 



V2 ' V2 
Imposing energy-momentum conservation at the partonic level via S'^{p + pA 

= k^ + , 
p\ = k++l+, 
PA,t = kt + k- 



1 

7^ 



1 

7f 



■t e + — e 



qt 

— f 



(A7) 
(A8) 



where the final hadron transverse momentum and rapidity are denoted by qt and rjh, and we used Zf = qt/lt- Note 
that light-cone momentum fraction Xg appears in the dipole forward scattering amplitude Np{bt, rt,Xg) whereas Xq is 
the fraction of the projectile proton (deuteron) carried by the incident quark, see Eq. (IA2[) . To derive an expression 
for the lower limit of integration in Eq. we note that < Xq < 1. Using the relation between the minus 
components of the four momenta given above, we get 
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The minimum value of occurs when Xq is maximum, i.e., Xq = 1. We then have 

9" 



which can be written in terms of the transverse momenta and rapidities of the final state hadron and photon as 



qt e"^ 



(All) 



We now consider the kinematics of single inclusive photon production cross section. The cross section is obtained 
from Eq. ([T|) by integrating over the final state quark momenta. This requires some care as we have now explicitly 
separated direct and fragmentation photons in Eq. ([TU)) . Again using Eqs. (IA21 IA31 IA4[) . we obtain the following 
relation, 

xq^^ (h e-"- + J^) , (A12) 



where opposite to Eq. (jASp . we avoided to introduce rjh and Zf. One can use the energy-momentum delta functions 
in Eq. (IA4[) to obtain the following relation 



I = Xq^/S/2- kt/V2e-'^-'. (A13) 
Now using the above relation and Eq. (|A12p . we obtain 

1 



Xq 



Xq S 



, ^? (A14) 



z 1 — z 

where the parameter z is the fraction of energy of parton carried away by photon and it is defined as follows 

P~ XqVS 



(A15) 



In case of direct photons with transverse momentum kt , one should shift momentum It It ^ kt in Eq. ()A14p (this 
is how we obtained the expersion Eq. (ITUl) '). Assuming that k « kt, we get 

where we have now used Xg to denote the light cone momentum fraction of the target carried by gluons for production 
of direct photons so as to distinguish it from the momentum fraction involved in production of fragmentation photons. 
In the later case, the integration variable k has been shifted twice. Implementing the shifts in Eq. (|A14I) and noting 
that the It integration in the fragmentation photon production cross section is dominated by its singularity at — > 
we get, for fragmentation photons, 

= (A17) 

Z Xqi^ 
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